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ABSTRACT 

In hierarchical models, where spheroidal galaxies are primarily produced via a 
continuous merging of disk galaxies, the number of intrinsically red systems at faint 
limits will be substantially lower than in "traditional" models where the bulk of star 
formation was completed at high redshifts. In this paper we analyse the optical-near- 
infrared colour distribution of a large flux-limited sample of field spheroidal galaxies 
identified morphologically from archival Hubble Space Telescope data. The Isi4 — HK' 
colour distribution for a sample jointly limited at /814 <23 mag and HK' <19.5 
mag is used to constrain their star formation history. We compare visual and au- 
tomated methods for selecting spheroidals from our deep HST images and, in both 
cases, detect a significant deficit of intrinsically red spheroidals relative to the predic- 
tions of high-redshift monolithic collapse models. However the overall space density of 
spheroidals (irrespective of colour) is not substantially different from that seen locally. 
Spectral synthesis modelling of our results suggests that high redshift spheroidals are 
dominated by evolved stellar populations polluted by some amount of subsidiary star 
formation. Despite its effect on the optical-infrared colour, this star formation prob- 
ably makes only a modest contribution to the overall stellar mass. We briefly discuss 
the implications of our results in the context of earlier predictions based on models 
where spheroidals assemble hierarchically. 
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1 INTRODUCTION 

The age distribution of elliptical galaxies is a controversial 
issue central to testing hierarchical models of galaxy forma- 
tion. The traditional viewpoint (Baade 1957, Bandage 1986) 
interprets the low specific angular momentum and high cen- 
tral densities of elliptical galaxies with their dissipationless 
formation at high redshift. In support of this viewpoint, ob- 
servers have cited the small scatter in the colour-magnitude 
relation for cluster spheroidals at low redshifts (Sandage & 
Visvanathan 1978, Bower et al 1992) and, more recently, 
such studies have been extended via HST imaging to high 
redshift clusters (Ellis et al 1997, Stanford et al 1997). Ex- 
amples of individual massive galaxies with established stel- 
lar populations can be found at quite significant redshifts 
(Dunlop 1997). 

In contrast, hierarchical models for the evolution of 
galaxies (Kauffmann et al 1996, Baugh et al 1996) pre- 



dict a late redshift of formation for most galactic-size ob- 
jects because of the need for gas cooling after the slow 
merger of dark matter halos. These models propose that 
most spheroidal galaxies are produced by subsequent merg- 
ers of these systems, the most massive examples of which 
accumulate since z :^1. Although examples of apparently 
old ellipticals can be found in clusters to quite high redshift, 
this may not be at variance with expectations for hierarchi- 
cal cold dark matter (CDM) models since clusters represent 
regions of high density where evolution might be accelerated 
(Governato et al 1998). By restricting evolutionary studies 
to high density regions, a high mean redshift of star forma- 
tion and homogeneous rest-frame UV colours would result; 
such characteristics would not be shared by the field popu- 
lation. 

Constraints on the evolution of field spheroidals derived 
from optical number counts as a function of morphology 
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(Glazebrook et al 1995, Im et al 1996, Abraham et al 1996a) 
are fairly weak, because of uncertainties in the local luminos- 
ity function. Nonetheless, there is growing evidence of differ- 
ential evolution when their properties are compared to their 
clustered counterparts. Using a modest field sample, Schade 
et al (1998) find a rest-frame scatter of S{U - V)^0.27 for 
distant bulge-dominated objects in the HST imaging sur- 
vey of CFRS/LDSS galaxies, which is significantly larger 
than the value of ^^0.07-0.10 found in cluster spheroidals 
at z :^0.55 by Ellis et al 1997. Likewise, in their study of 
a small sample of galaxies of known redshift in the Hubble 
Deep Field (HDF), Abraham et al (1998) found a significant 
fraction (~40%) of distant ellipticals showed a dispersion in 
their internal colours indicating they had suffered recent star 
formation possibly arising from dynamical perturbations. 

Less direct evidence for evolution in the field spheroidal 
population has been claimed from observations which at- 
tempt to isolate early-type systems based on predicted 
colours, rather than morphology. Kauffmann et al (1995) 
claimed evidence for a strong drop in the volume density of 
early- type galaxies via a V/Vmax analysis of colour-selected 
galaxies in the Canada- France Redshift Survey (CFRS) sam- 
ple (Lilly et al 1995). Their claim remains controversial 
(Totani & Yoshii 1998, Im & Castertano 1998) because of 
the difficulty of isolating a robust sample of field spheroidals 
from V — I colour alone (c.f. Schade et al 1998), and the 
discrepancies noted between their analyses and those con- 
ducted by the CFRS team. 

In addition to small sample sizes, a weakness in most 
studies of high redshift spheroidals has been the paucity of 
infrared data. As shown by numerous authors (eg. Chariot & 
Silk 1994), near-IR observations are crucial for understand- 
ing the star formation history of distant galaxies, because at 
high redshifts optical data can be severely affected by both 
dust and relatively minor episodes of star- format ion. Recog- 
nizing these deficiencies, Moustakas et al (1997) and Glaze- 
brook et al (1998) have studied the optica- infrared colours of 
small samples of of morphologically-selected galaxies. Zepf 
(1997) and Barger et al (1998) discussed the extent of the 
red tail in the optical-IR colour distribution of HDF galax- 
ies. Defining this tail (Veoe-K >7 and Iqia-K >4) in the 
context of evolutionary tracks defined by Bruzual & Char- 
lot's (1993) evolutionary models, they found few sources in 
areas of multicolour space corresponding to high redshift 
passively-evolving spheroidals. 

The ultimate verification of a continued production of 
field ellipticals as required in hierarchical models would be 
the observation of a decrease with redshift in their comov- 
ing space density. Such a test requires a large sample of 
morphologically-selected ellipticals from which the luminos- 
ity function can be constructed as a function of redshift. By 
probing faint in a few deep fields, Zepf (1997) and Barger 
et al (1998) were unable to take advantage of the source 
morphology; constraints derived from these surveys relate to 
the entire population. Moreover, there is little hope in the 
immediate term of securing spectroscopic redshifts for such 
faint samples. The alternative adopted here is to combine 
shallower near-infrared imaging with more extensive HST 
archival imaging data, allowing us to isolate a larger sample 
of brighter, morphologically- selected spheroidals where, ulti- 
mately, redshifts and spectroscopic diagnostics will become 
possible. Our interim objective here is to analyse the optical- 



infrared colour distribution of faint spheroidals which we will 
demonstrate already provides valuable constraints on a pos- 
sible early epoch of star formation. 

A plan of the paper follows. In §2.1 we discuss the avail- 
able HST data and review procedures for selecting morpho- 
logical spheroidals from the images. In §2.2 we discuss the 
corresponding ground-based infrared imaging programme 
and the reduction of that data. The merging of these data 
to form the final catalogue is described in §2.3. In §3 we dis- 
cuss the optical-infrared colour distribution for our sample 
in the context of predictions based on simple star forma- 
tion histories and consider the redshift distribution of our 
sample for which limited data is available. We also exam- 
ine constraints based on deeper data available within the 
Hubble Deep Field. In §4 we summarise our conclusions. 



2 CONSTRUCTION OF THE CATALOGUE 
2.1 THE HST SAMPLE 

In searching the HST archive for suitable fields, we adopted 
a minimum / F814W-band exposure time of 2500 sec and 
a Galactic latitude of |6| = 19° so that stellar contamination 
would not be a major concern. These criteria led to 48 fields 
accessible from the Mauna Kea Observatory comprising a to- 
tal area of 0.0625 deg^(225 arcmin^). Table 1 lists the fields 
adopted, including several for which limited redshift data 
is available e.g. the HDF and its flanking fields (Williams 
et al 1996), the Groth strip (Groth et al 1994) and the 
CFRS/LDSS survey fields (Brinchmann et al 1997). F606W 
imaging is available for 25 of the fields in Table 1. Object 
selection and photometry for each field was performed using 
the SExtractor package (Bertin & Arnouts 1996). Although 
the detection limit varies from field to field, the Iqia -band 
data is always complete to ^ 24 mag and the Veoe -band to 
- 25 mag0. 

The morphologies of galaxies in the sample were inves- 
tigated independently using visual classifications made by 
one of us (RSE), and automated classifications based on the 
central concentration (C) and asymmetry [A) parameters 
defined in Abraham et al. (1996b). In the case of visual clas- 
sifications we adoped the MDS scheme defining spheroidal 
to include E: E/SO: SO and SO/a (MDS types 0,1,2). As 
shown below, the visual and automated classifications com- 
pare quite favourably, with particularly satisfactory agree- 
ment for the regular spheroidal galaxies that are the focus 
of this paper. 

The appropriate limiting magnitude of our survey is 
set by that at which we believe we can robustly isolate 
spheroidal galaxies from compact HII galaxies, stars and 
bulge-dominated spirals. The Medium Deep Survey (MDS) 
analyses adopted a limiting magnitude for morphological 
classification (using nine classification bins) of Iqia =22 mag, 
although some MDS papers extended this further to Iqia 
=23 mag (see Windhorst et al 1996 for a summary). In Abra- 
ham et al (1996b) and Brinchmann et al (1997), HST data 
similar to that in the present paper was also used to classify 

^ These and subsequent detection limits refer to near-total mag- 
nitudes in the Vega system based on profile fitting within the 
SExtractor package ('mfoest')- 
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Figure 1. A selection of visually- classified spheroidal galaxies sorted by /8i4 magnitude, selected from the HST archive. Each panel 
represents 10 arcsec on a side. 



galaxies to /8i4 =22 mag. However, by restricting our anal- 
ysis in the present paper to spheroidal systems, we are able 
to extend classifications to slightly deeper limits (/8i4 =23.0 
mag). This is possible because the chief diagnostic for dis- 
criminating spheroidals is central concentration, rather than 
asymmetry which is sensitive to lower surface brightness fea- 
tures. Because the classifications based on A and C are ob- 
jective, the classification limits for the present dataset have 
been investigated using simulations, as described below. 
Figure 1 shows a typical set of morphologically- 



identified spheroidals at various magnitudes down to Iqia 
=23 mag. Figure 2 compares the A — C and visual mor- 
phological distributions at a range of magnitude intervals, 
down to the limits of our survey. The demarcation between 
early and late-types on the basis of A and C is made us- 
ing bright galaxies (/8i4 < 20 mag) and shifted slightly as 
a function of magnitude on the basis of simulations made 
using the IRAF package art data, which model the apparent 
change in the central concentration of an r^^^ law elliptical 
galaxy as a function of decreasing signal-to-noise. Random 



© 0000 RAS, MNRAS 000, 0-0 



4 F. Menanteau et al 



-0.5 



o 



-1.5 



I8i4<20 



o ^ o 
o 

Of 



o 
o 



-1 -0. 



-0.6 -0.4 -0.2 
Log(C) 



3- 

o 



-0.5 



-1.5 - 




-1 -0. 



-0.6 -0.4 -0.2 
Log(C) 



-0.5 



o 



■1.5 



-2 



21<I,^,<22 



o 



O (Q ^ 0\3 

O o-^o* 



O •qOO o 




-0.6 -0.4 
Log(C) 



-0.2 



-0.5 



o 



1.5 



-2 



22<Iaj4<23 



n CP 



o 



o 
o 

o o 







_l ^ 



-O.f 



-0.6 -0.4 
Log(C) 



-0.2 



Figure 2. T/ie distribution of asymmetry and concentration for visually- classified galaxies a) [top left] /8i4 < 20 mag^ h) [top right] 
20<IsiA <21 mag, c) [bottom left] 21<l8i4 <22 mag (the limit achieved with the Medium Deep Survey) and (d) [bottom right] 22<Isia 
<23 mag. The dashed lines represent optimal boundaries for the separation of spheroidal galaxies (solid circles) from spirals and 
irregulars (open circles). Compact objects are shown as crosses. The reference boundary line is defined using the visually identified 
sample at IgiA <20 and then shifted slightly as a function of magnitude on the basis of simulations. The expected la RMS errors on 
measures of central concentration are shown by the error bar on the lower right portion of each panel. 



errors on central concentration are also determined on the 
basis of simulations, and representative error bars are shown 
in Figure 2. The general agreement between the visual and 
automated classifications is remarkably good, particularly 
to /814 =22 mag. Between IgiA =22 and 23 mag the agree- 
ment worsens, mostly because of the great increase in the 
number of visually- classified "compact" systems. We define 
compacts to be those systems where there is no clear distinc- 
tion between small early- type galaxies, faint stars and/or 
HII regions. 

In order to quantify the concordance between the vi- 
sual and automated classifications, the A — C distribution 
was analysed using a statistical bootstrap technique (Efron 
& Tibshirani 1993). The A — C distribution was resampled 
500 times in order to determine the uncertainties in both the 
number of systems classified as early- type, and the uncer- 
tainties on the colour distribution for these systems. These 
measurements will be discussed further below in §3.3. 

The somewhat larger number (323 vs 266) of A/C- 



classified early-types relative to the visually classified galax- 
ies is significant at the 3a level. However, if the compact 
systems are included in the tally of visually classified early- 
type systems, then the number of visually and A/C classified 
ellipticals agree closely (to within 1 sigma). It is clear that 
the distinction between compact galaxies and early-type sys- 
tems is an important consideration when determining the 
number counts of early type systems at the faint limits of 
these data. However, it is perhaps worth noting at this stage 
that another bootstrap analysis (presented in §3.3) shows 
that the uncertainty introduced by compact systems into 
the number counts at 22 < / < 23 does not manifest it- 
self as a large uncertainty in the colour histograms of the 
early- type population. 



2.2 GROUND-BASED INFRARED IMAGING 

Although some of the fields in Table 1 have Veoe and IgiA 
HST data, such a wavelength baseline is not very useful 
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Figure 3. Transmission curve for the notched HK' filter (solid 
line) compared with that for the standard H (dashed line) and 
K' (long-dashed line) passbands 

at thesed depths. As discussed by Moustakas et al (1997) 
and Zepf (1997), the addition of infrared photometry is es- 
pecially helpful in distinguishing between passively-evolving 
systems and those undergoing active star formation, regard- 
less of redshift, primarily because of its reduced sensitivity 
to K-dimming, small amounts of star formation and dust 
reddening. 

Our infrared imaging was mainly conducted using the 
QUIRC 1024^ infrared imager on the University of Hawaii 
2.2-m telescope. The log of observations is summarised in 
Table 2. In order to improve the observing efficiency in secur- 
ing deep infrared photometry for a large number of WFPC-2 
fields, we used the notched H + K' l.S/xm filter (which we 
refer to hereafter as the HK' filter) (Wainscoat & Cowie 
1998, Figure 3) which offers a gain in sensitivity of typically 
a factor of ^2 over a conventional K' filter. At the f/10 focus 
of the UH 2.2m, the field of view is 193'' x 193'' with a scale 
of 0.1886" pixel" ^ ensuring that the 3 WFPC2 chips can be 
comfortably contained within a single exposure. 

Each HK' exposure was composed of 13 sub-exposures 
of ~100 sec duration (depending on the background level) 
spatially-shifted by increments of 5-20 arcsec in all direc- 
tions. This dithering pattern was repeated 2-3 times dur- 
ing the exposure. The data was processed using median sky 
images generated from the disregistered exposures and cal- 
ibrated using the UKIRT faint standards system. Most of 
the data was taken under photometric conditions; deep non- 
photometric data was calibrated via repeated short expo- 
sures taken in good conditions. The limiting magnitude of 
the infrared data varies slightly from field to field and is 
deepest for the HDF and flanking fields which were taken in 
a separate campaign (Barger et al 1998). 

In order to determine the detection limit of our HK' 
data we performed extensive Monte Carlo simulations. Us- 
ing the IRAF art data package we created simulated data 
sets, which were subsequently analysed using the same ex- 
traction and measurement methods as for the real data. 
With the task mkobjects we generated artificial galaxies as- 
suming exponential disk profiles with no internal absorption 
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Figure 4. Statistical completeness of the UH 2.2m QUIRC data 
for the EST archival fields and the HDF as determined by pro- 
cedures described in the text. 

for spirals and de Vaucouleurs profiles for spheroidals. The 
profile scales were chosen to be magnitude-dependent con- 
verging to the image seeing at faint limits. Figure 4 shows 
the results of this exercise. The 80% completeness limit for 
spheroidals is HK' =19.5 mag for most of the survey extend- 
ing to HK'=20.b mag for the HDF and flanking fields. 

2.3 COMPLETENESS OF THE COMBINED 
OPTICAL-INFRARED CATALOGUE 

The final photometric catalogue of spheroidals was obtained 
by matching the HST /8i4-band and the ground-based IR 
SExtractor catalogues using the adopted magnitude limits 
of HK' <19.5 mag and /8i4 <23.0 mag. In the final matched 
catalogue, we retained the SExtractor ^rrihest magnitudes 
but measured / — HK' colours within a fixed 3 arcsec diam- 
eter. This aperture size, together with the fairly good seeing 
of the IR data, ensures that when calculating colours we 
are looking at the same physical region of the galaxy. Of the 
818 sources in the matched catalogue, 266 systems were visu- 
ally classified as spheroidals (defined to be one of 'E, E/SO, 
SO, or so/a' in the MDS scheme) and 50 as compact ob- 
jects. Automated classifications result in 323 sources classed 
as spheroidals (with no distinction between spheroidals and 
compacts). 

Clearly the joint selection by Iqia and HK' necessary 
to exploit HST's morphological capabilities and establish 
optical-infrared colours could lead to complications when 
interpreting / — HK' colour distributions. As a major moti- 
vation for this study is to identify as completely as possible 
the extent of any red tail in the colour distribution, incom- 
pleteness caused by the various magnitude limits is an im- 
portant concern. Figure 5 shows that, within the optical, 
morphologically- selected sample with /8i4 <23 mag, virtu- 
ally all of the HK' <19.5 mag sample can be matched; only 
a small fraction (18/818=2.2%) of red / - HK' >3.5-5 mag 
objects are missed. We return to the nature of these sources 
in §3.3. 
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Figure 5. Colour-magnitude distribution for the visually- 
classified catalogue limited at HK' =19.5 mag and /8i4 =23 mag. 
Solid points correspond to spheroidals, crosses to compacts and 
circles to the remaining spirals and late type galaxies. The joint 
selection in /8i4 and HK' implies a small fraction (<3%) of 
HK' -selected objects are not contained within the HST sample. 
These objects are shown as stars and arrows as lower limit when 
no detection was possible 



3 ANALYSIS 
3.1 Strategy 

Our analysis is motivated by the two principal differences we 
might expect between models where ellipticals underwent a 
strong initial burst of activity with subsequent passive evo- 
lution (which we will term the 'monolithic collapse' model) 
and those associated with a hierarchical assembly of ellipti- 
cals from the dynamical merger of gas-rich disks (Baugh et 
al 1996). We recognise at the outset that these models repre- 
sent extreme alternatives with a continuum of intermediate 
possibilities (c.f. Peacock et al 1998; Jimenez et al. 1998). 
Our strategy in this paper, however, will be to discuss our 
field elliptical data in the context of the simplest models pro- 
posed to explain the star formation history of distant cluster 
ellipticals (Ellis et al 1997, van Dokkum et al 1998). More 
elaborate analyses are reserved until spectroscopic data is 
available for a large sample. 

Firstly, in the monolithic collapse model, the comoving 
number density of luminous ellipticals should be conserved 
to the formation redshift (say, z ^3-5), whereas in hierar- 
chical models we can expect some decline in number density 
at moderate redshift depending on the cosmological model 
and other structure formation parameters (Kauffmann et 
al 1996, Kauffmann & Chariot 1998a). Such a change in 
the absolute number density would be difficult to convinc- 
ingly detect without spectroscopic data. The number of faint 
HST- identified ellipticals has been discussed by Glazebrook 
et al (1995), Driver et al (1995), Abraham et al (1996b) 
and Im et al (1996) with fairly inconclusive results because 
of uncertainties arising from those in the local luminosity 
function used to make predictions (Marzke et al 1998). 

Secondly, there will be a redshift-dependent colour shift 



Figure 6. Colour histograms for various morphological samples 
within our adopted limits of HK' <19.5 and /8i4 =23. The solid 
line refers to the visually -classified spheroidals (MDS types 0,1,2) 
and the dashed- dot line to those spheroidals identified on the basis 
of their asymmetry and concentration indices. The long dashed 
line refers to visually- classed compacts (type=-l) and the short 
dashed line refers to the remainder (types >3). 

associated with merger-driven star formation in the hierar- 
chical models whereas, for the monolithic case, the sources 
will follow the passive evolution prediction. Kauffmann et al 
(1996) initially claimed that both signatures would combine 
in the hierarchical picture to produce a factor 3 reduction 
in the abundance of passively-evolving sources by z ^1, but 
a more recent analysis (Kauffmann & Chariot 1998b) shows 
that the decline is dependent on the input parameters. For 
example in a model with non-zero cosmological constant (the 
so-called 'ACDM'), little decline is expected until beyond 
z:^l. 

In contemplating these hypotheses in the context of our 
HST data, it must be remembered that although HST can be 
used very effectively to isolate spheroidals morphologically 
to /814 =23 mag (representing a considerable advance on 
earlier colour-selected ground-based samples which could be 
contaminated by dusty later types), in the case of merger 
models, the predictions will depend critically on the time 
taken before a system becomes a recognisable spheroidal. 
However, any hypothesis which postulates a constant co- 
moving number density of well-established spheroidals is 
well suited for comparison with our data, the outcome be- 
ing important constraints on the past star formation history 
and luminosity evolution. 

3.2 Colour Distributions 

Figure 6 shows /8i4 -HK' colour histograms for both the 
visual and A/C-selected spheroidals alongside those for the 
compacts and the remainder. The automated and visual cat- 
alogues have nearly identical colour distributions, confirming 
earlier tests on the reliability of the automated classifier. In 
fact, the differences between the automated and visually de- 
fined histograms are almost completely attributable to the 
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compact systems, which cannot be segregated from other 
early- types on the basis of central concentration. The colour 
histogram for compacts spans the range seen for early-type 
galaxies, with a peak slightly red ward of that for visually 
classified early- types. It is clear from the similarity between 
the colour histograms for visual and automated classifica- 
tions that contamination of spheroidals by compacts (ex- 
pected in the automated catalogue) does not pose a signifi- 
cant uncertainty in determining the colour distribution. 

The histogram of colours for late-type galaxies peaks at 
nearly the same colour as that for the early- types, which at 
first seems somewhat surprising. As we will later see, this 
is largely a refiection of the very wide redshift range in- 
volved. However, the distribution for spirals and later types 
is skewed toward the blue, although redward of I—HK' =2.5 
mag the shapes of the distributions are similar (see also §3.5 
below) . 

3.3 Single Burst Model Predictions 

Figure 7 compares the observed colour histograms with a 
range of model predictions based on the GISSEL96 spectral 
synthesis code (Bruzual & Chariot 1996) for a range of star- 
formation histories. Observed and predicted total counts for 
each of the models are also given in Table 3. At this stage we 
concentrate on 'single burst' or 'monolithic collapse' models 
which conserve the comoving number density at all epochs, 
and defer discussion of alternative scenarios until §3.5. 

Our model predictions take into account the joint Iqia 
and HK' selection criteria for our sample, and are based 
on the present-day optical E/SO luminosity function (LF) 
from Pozzetti et al. (1996), ie. a standard Schechter function 
with = 0.95 X 10"^ Mpc"^, M^^ = -20.87 and a faint-end 
slope of a = —0.48. When making model predictions, this lu- 
minosity function is tranformed into one appropriate for the 
/814 photometric band via a single colour shift, resulting in 
^^814 ~ —23.12. For comparison, we also show predictions 
assuming a suitably transformed luminosity function with a 
fiat faint-end slope (a = -1) and = 0.55 x 10~^Mpc~^ 
as suggested by Marzke et at (1998). Throughout this paper 
we adopt Hq — 50 Km s~^Mpc~^. Given the elementary 
nature of the comparisons currently possible, and the fact 
that the expected dispersion in / — HK' from the present- 
day colour- luminosity relation is minimal, we have avoided 
the temptation to model a distribution of metallicities within 
the galaxy population, preferring instead to explore the ef- 
fects of fixing the metallicity of the entire population to a 
single value within a large range (40%-250% solar) in the 
simple predictions discussed below. Other variables in the 
single burst hypothesis include the redshift of formation, zf 
(fixed at Zf — 5), the burst- duration (represented as a top 
hat function of width 1.0 Gyr) and the cosmological parame- 
ters (Qm and Q.a)- As shown in Appendix A, the luminosity 
weighted metallicities of the present sample are not biased 
strongly by the limiting isophotes of the our observations, 
and fair comparisons can be made using individual single- 
met allicity tracks over a broad range of redshifts. 

Clearly the most important input parameter in the 
model predictions shown in Figure 7 (summarized in Ta- 
ble 3) is the metallicity. Although our spheroidals are almost 
exclusively luminous (> L*) galaxies which, in the context 
of single-burst models would imply a metallicity of at least 



solar (cf. Appendix A, and Arimoto et al 1997), here we will 
explore the possibility that part of the colour distribution 
could arise from a wider metallicity range than that found 
locally. 

3.3.1 A Deficit of Red Spheroidals 

The predicted colour distributions show a characteristic 
asymmetry. This is caused by the blending of the Iqia—HK' 
K- correction and the passive bluing of systems to z ^1.5. 
In contrast, the observations reveal a clear excess of blue 
(/814 —HK' <2 mag) spheroidals not predicted by even 
the lowest metallicity models. More significantly, solar and 
super-solar metallicity models over-predict the extent of the 
red tail in the colour distribution. Both discrepancies are 
consistent with recent star- format ion in our sample of faint 
spheroidals. In order to quantify these discrepancies, we per- 
formed a Kolmogorov-Smirnov test (K-S) to check whether 
the observed spheroidal colours could be drawn from any 
of the model distributions (allowing a measurement error 
<^i-HK' = 0.25 mag). In all cases the observed distribution 
differs from the models at a confidence level higher than 
99.99%. Evidently monolithic collapse models with constant 
co-moving density fail to reproduce the colour distribution 
of high redshift spheroidals. 

It will be convenient in the following to quantify the 
strength of the red tail in the colour distribution by defin- 
ing a "red fraction excess", shown in Table 3, constructed 
as the ratio of the predicted number of early type systems 
with /814 — HK' > 3.0 mag to the observed number. The 
statistical uncertainties on the red fraction excess in this ta- 
ble are based on 500 bootstrap resamplings of the original 
catalogue, each realization of which was subjected to the 
same selection criteria applied to the original data. 

As discussed by many authors (Glazebrook et al 1995, 
Marzke et al 1998), the absolute numbers depends sensi- 
tively on the normalisation and shape of the local LF. Table 
3 includes a summary for the range in LF parameters dis- 
cussed earlier. For a declining faint-end slope (a = —0.48) 
and solar- and- above metallicity, the red fraction excess is 
more than 5 times that observed. Adopting a metallicity 
substantially below solar results in closer agreement but as- 
suming such metallicities for the entire population may be 
unreasonable given local values (see Appendix, and Arimoto 
et al 1997). Even so, the red fraction excess is only reduced 
from 8 to ~ 3 if the adopted metallicity is varied between 
250% solar and 40% solar. Models with a fiat faint-end slope 
(a = —1) improve the agreement further and in the very 
lowest metallicity model with Qm = there is almost no 
deficit. 

In Table 3 we have also included models with Qa > 
for both slopes of the LF. The effect of Qa is to increase 
the apparent deficit of red spheroidals (because of the rapid 
increase in the differential volume element with redshift for 
Qa > cosmologies at z < 1). 

Although models where the redshift of the initial burst, 
Zf, is reduced to 2; = 3 are not shown in the table, these 
do not result in significant changes in the above discussions. 
We conclude that we cannot reconcile the number of galax- 
ies in the red end of the observed colour histogram to the 
corresponding predictions of a constant co-moving density 
high-redshift monolithic collapse model. Alternative scenar- 
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Figure 7. The colour distribution of visually- classified spheroidals with various single burst models (see legend) compared to the observed 
number, represented by solid histogram (see text for further assumptions). 



ios, which may explain the relatively blue colours of some 
observed spheroidals, will be considered in §3.5. 



3.3.2 A Declining Number of High Redshift Spheroidals? 

While monolithic collapse models fail to reproduce the ob- 
served colour distributions, Table 3 indicates that the over- 
all number is in reasonable agreement. Specifically, for a low 
Qm and ^^a^O, the Marzke et al. luminosity function and 
luminosity- weighted metallicities of solar and above, we see 
no significant evolution in the space density of spheroidals. 
For the currently popular spatially-flat universe with low 
Qm and high Qa (Perlmutter et al 1999), the data imply 
a deficit of no more than 30%. Stronger evolution 60% 
decline) would occur if we adopted the Pozetti et al. lumi- 
nosity function. We therefore conclude that the colour offset 
described earlier is more likely the result of star- format ion 
activity in well-formed spheroidals at high redshifts rather 
than evidence for evolution in their space density. 

At this point, we return to the nature of those 18 sources 
identified in the infrared images which are fainter than /8i4 
=23 mag. Although they could formally be included in the 
colour distributions, they are too faint in the WFPC-2 im- 
ages for reliable morphological classification. A montage of 
these sources is shown in Figure 8. At most 3 of the 18 are 



possible spheroidals. As such, their addition to the colour 
distribution would have a negligible impact on conclusions 
drawn from Figure 7. 



3.4 Constraints from Redshift Distributions 

While our principal conclusions are based on the enlarged 
size of our HST sample and addition of infrared photometry 
compared to ear Her work, it is interesting to consider what 
can be learnt from the (incomplete) redshift data currently 
available for our sample. We have collated the published 
spectroscopic redshift data from the CFRS/LDSS surveys 
(Brinchmann et al 1998), the MDS survey (Glazebrook et 
al 1998) and the HDF and its fianking fields (tabulated by 
Cowie 1997) and matched these with our HK' -selected sam- 
ple. In total, 97 of our galaxies have published redshifts. As 
the bulk of these surveys were themselves magnitude- limited 
in /, the magnitude and colour distribution of this subset 
should be representative of that for an appropriate subset 
of our primary photometric sample. 

Figure 9 (upper panel) shows the colour-redshift rela- 
tion for the 97 objects in the subsample with redshift in- 
formation. Also shown are evolutionary predictions based 
on spectral synthesis models adopting ranges in metallicity 
and star- format ion history as before. The 47 spheroidals in 
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Figure 8. A montage of WFPC-2 images for those red sources contained within the HK' <19.5 sample hut for which IsiA. >23. 



this set are clearly redder at a given redshift than their spiral 
and later-type counterparts and span a wide redshift range 
with median value 'z ^0.7. However, as Schade et al (1998) 
discussed in the context of Veoe -/8i4 colours for their smaller 
sample of HST field galaxies, there is some overlap between 
the classes at a given redshift. The colour scatter for morpho- 
logical spheroidals appears somewhat larger than the ^0.2 
mag dispersion expected from slope of the infrared-optical 
colour-luminosity relations for early-type systems^. 

In the specific case of the HDF, it is interesting to ex- 
ploit the increased depth of both the HK' data and the HST 
optical morphologies (Abraham et al 1996a) as well as to 
consider the abundant photometric redshift estimates. For 
this purpose we constructed a 19. 5< HK' <20.5 mag exten- 
sion to our HDF sample, with morphological classifications 
from the deep (/8i4 < 25 mag) morphological catalogue of 
van den Bergh (1996). For this sample we can take advan- 
tage of the apparently rather good precision in photomet- 
ric redshift estimates for early- type galaxies (Connolly et al 
1997, Wang et al 1998)|. By going deeper in the HDF we 
extend our sample by 20 objects, of which 8 are classed visu- 
ally as E/SOs (none are compact). Adding this extension to 
those HDF galaxies already in our catalog, the correspond- 
ing numbers with HK' < 20.5 mag; /8i4 < 25 mag become 
50 and 26 respectively. The colour-redshift relation for the 
combined HDF sample is shown in Figure 9 (lower panel). 
For the same sample, Figure 10 shows the colour-magnitude 
diagram of the visually classified E/SOs. The inset shows 
the colour histogram and the arrow indicates the peak in 
the distribution for the primary sample. 



t Note however that the slope of the infrared-optical colour- 
magnitude relation for early-type systems is rather uncertain, par- 
ticularly for I—K. On the basis of quite strongly model-dependent 
conversions based on the cluster V—K data of Bower et al. (1992), 
Peletier &; de Grijs (1998) obtain a slope of -0.0438 ± 0.0041 for 
the I — K slope of local early- type systems, using the models 
presented in Vazdekis et al. (1996). 

t We note that the good accuracy in photometric redshifts for 
these galaxies appears to be due to the presence of strong con- 
tinuum features which are well-explored with the addition of the 
Kitt Peak JHK photometry to the HDF filter bands. 



As expected, the peak of the HDF colour histogram in 
Figure 10 lies redward of the colour histogram for our entire 
sample (by ~ 0.2 mag). But the redshift data in Figure 9 
makes it clear that this peak is still substantially bluer than 
expected for the simple monolithic collapse model at solar 
metallicity. For the 26 HDF spheroidals, spectroscopic red- 
shifts are available for 19, the rest being photometric. Im- 
portantly, the spectroscopically-confirmed galaxies include 3 
ellipticals beyond z:^0.9, all of which are substantially bluer 
than the passive evolution predictions. While based on small 
numbers of galaxies, the figure lends strong support to the 
conclusions of the previous subsection, particularly when it 
is realised there is an in-built bias in favour of photomet- 
ric redshifts matching the passively-evolving spectral energy 
distributions. 

These conclusions based on the HDF are consistent with 
those of Zepf (1997) and Barger et al (1998) who analysed 
optical- infrared colours of much fainter sources without tak- 
ing into account morphological and redshift information. 



3.5 Alternative Star Formation Histories 

The single burst models ruled out by the colours of sper- 
oidals in the previous sections are idealised representations 
of spheroidal history. We now consider alternative histories 
which could be more consistent with our various datasets. 

At its most fundamental level, the deficit of red 
spheroidals at faint limits appears to eliminate models with 
very short epochs of star formation at high redshifts fol- 
lowed by long quiescent periods. In the context of modelling 
distant red radio galaxies, Peacock et al (1998) have shows 
that models with continuous star formation truncated at 
later times avoid the peak luminosities associated with ini- 
tial bursts and can produce a significant bluing at redshifts 
where the red tail would otherwise be seen^. In an important 



^ On the basis of model predictions used to calculate the density 
of post-starburst strong" systems seen in clusters, Couch & 
Sharpies (1986), Barger et al. (1996) and Abraham et al. (1996c) 
note that a sharp truncation in the star-formation rate of ac- 
tively star-forming systems results in a synchronization of optical 
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Figure 10. The colour- magnitude diagram of visually- classed objects within the deeper sample possible solely for the HDF (IsiA < 23^ 
HK' <20.5). The inset histogram shows the IgiA — HK' colour distribution for visually- classed spheroidal within the deeper sample. 
The arrow shows the peak of the distribution in the primary sample. Symbols are as for Figure 9 



recent paper, Jimenez et al. (1998) argue that monolithic 
models have been rejected prematurely by some authors: 
only extreme scenarios with very short duration bursts (eg. 
10^ Myr) of star- format ion followed by absolute quiescence 
can be ruled out, while bursts with fairly low- levels of ex- 
tended star- format ion activity may still be compatible with 
the observed data. On the basis of a simple one-dimension 
chemo- dynamical model for the evolution of spheroids, these 
authors predict that the integrated star- format ion history of 
early-type galaxies should resemble a quasi-monolithic col- 
lapse model. A concrete prediction of this model is that the 
bulk of the star- format ion in high-z spheroids is occuring 
near their cores, an effect that appears to have been seen 
(Abraham et al. 1999), and which may be responsible for 
the blue colours of some ellipticals in our sample. 

In the light of the above, it is therefore interesting to 
calculate the amount of recent star formation which must 
be added to a single burst model in order for the model 



colours with those expected of early type galaxies after only ~ 1.5 
Gyr, so it is not too surprising that truncated star formation his- 
tories and monolithic collapse models predict similar colours pro- 
vided the epoch of observation is several Gyr after the period of 
truncation. 



to match the observed /8i4 — HK' colours. We seek to de- 
termine the timescales over which such a system is bluer 
than the passive evolution model by at least Scoiour — 
(/sM - HK')modei - (hiA " HK')ohs ^ 0.3 mag, i.e. consis- 
tent with the typical colour offset from the solar metallicity 
tracks shown in Figure 9. To model this we added a short 
burst of star- format ion (of duration 0.1 Gyr, forming 15% 
of the stellar mass) occurring at z = Zhurst to an underly- 
ing solar metallicity monolithic collapse model at a given 
redshift using the population synthesis models of Bruzual & 
Chariot (1996). We then computed the redshift range over 
which the burst at Zhurst results in 6coior > 0.3 for a given 
redshift of observation, Zo- 

The result of this exercise is shown in Figure 11 for 
several redshifts of observation zq. The range in redshift 
space where a 15% burst produces a blueing 6coiour > 0.3 
is shown as the shaded area. (For example, a galaxy ob- 
served at redshift Zo = 1 would be at least 0.3 mag bluer 
than the monolithic collapse models between points A and 
B in the figure. At Zo = 0.5 a burst would have to occur at 
0.5 < Zhurst < 0.65, while for Zo — 0.7 it would have to oc- 
cur in the range Zhurst ^ 0.7 — 0.9. Since almost all points in 
Figure 9 lie blueward of the solar metallicity model tracks, 
and since galaxies in our sample lying in the redshift range 
< z < 1 have "memory" of bursts over 5z ^ 0.1 — 0.2, 
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Figure 11. The colour difference between a 1.0 Gyr single burst and the same model plus a short burst of star formation as a function 
of the redshift, z = Zfj^rst^ where the burst is added is show for a series of redshift of observation (zq = 0.2,0.5,0.7, 1.0, 1.5 and 2.0). 
The horizontal dashed line represent Scoiour = 0-3 and the shaded area shows the range in zijurst where 6coior ^ 0-3 for each Zq. 



it seems improbable that moderate intensity (ie. 15% mass) 
single burst events can explain the colour offsets shown in 
Figure 9. It seems more probable that the duty-cycle of star- 
formation is more extended, indicative of either a low level of 
roughly continuous star- format ion underlying the old stellar 
populations, or perhaps of a succession of lower mass bursts. 

The preceeding analysis indicates the extent to which a 
modest "polluting" star-forming population superposed on a 
dominant old stellar population reconciles our observations 
with traditional monolithic collapse scenarios. In contrast to 
this, it is interesting to consider what sort of hierarchical for- 
mation models may also be consistent with the present data. 
The rather mild density evolution in our sample is in sharp 
disagreement with the predicted factor of three decline in the 
abundance of spheroidals at z ^ 1, based on the present gen- 
eration of "semi-analytical" models in high-density, matter- 
dominated cosmologies (eg. Kauffmann & Chariot 1998a). 
However, in a flat A-CDM model (with Qm = 0.3, Qa = 0.7) 
the decUne in the abundance of spheroidals is only 30% at 
z ^ 1 (Kauffmann & Chariot 1998b), and may be consistent 
with the present observations. In the latter model the oft- 
quoted factor of three decline in the space density of ellipti- 
cals occurs at z = 2 instead of z = 1, so future work extend- 
ing the present sample to higher redshifts may allow us to 
distinguish between "extended" monolithic collapse scenar- 



ios and A-CDM hierarchical models. However, more detailed 
modelling in this particular case, e.g. of the colours distribu- 
tion, is beyond the scope of this paper. If star- format ion in 
extended monolithic-collapse scenarios is centrally concen- 
trated (as suggested by the simple one-dimensional models 
of Jimenez et al. 1998), the two scenarios may also be distin- 
guishable at lower redshifts on the basis of resolved colour 
distributions (Menanteau et al. 1999, in preparation). 



4 CONCLUSIONS 

We have constructed a new catalogue of ~300 faint field 
spheroidal galaxies using HST images for morphology to a 
limit of /814 =23 mag. Follow-up infrared photometry has 
enabled us to consider the optical-infrared colour distribu- 
tion of an infrared magnitude limited sample. We have mod- 
elled expected colours using various star formation histories, 
metallicities and cosmologies. For a limited subset, spectro- 
scopic redshift data is available and within the HDF it is 
possible to construct a deeper catalogue. 

Our main results can be summarised as follows: 
• There is little evidence for strong evolution in the 
space density of luminous field spheroidal systems out to z ^ 
1. Within the uncertainties introduced by poorly constrained 
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colours of high-redshift spheroidals in the present sample 
are consistent with the predictions of "extended" monolithic 
collapse scenarios, in which the existing star-formation pol- 
lutes the colours of a dominant, underlying old stellar pop- 
ulation. The data is not consistent with the predictions of 
semi-analytical hierarchical models in high-density, matter- 
dominated cosmologies. However, the observed weak density 
evolution may be consistent with the predictions of hierar- 
chical A-CDM models (Kauffmann & Chariot 1998b). 

Spectroscopic redshifts for a complete sub-sample of our 
catalogue will enable models for the star-formation history of 
high redshift spheroidals to be rigorously tested, eliminating 
some of the ambiguities present in the current analysis. 




Figure 9. The colour-redshift relation for that subset of the 
HK' <19.5 mag sample with published spectroscopic data (large 
symbols), and for a deeper HK' <20.5 mag sample within the 
HDF for which photometric redshift data (smaller symbols) is 
available (Wang et al 1998). Symbols refer to visual classifica- 
tions. The solid lines represents a passively evolving elliptical 
(single burst of duration 0.1 Gyr) at zp=5 for different metal- 
licities. Dashed lines show the corresponding models for zf=3. 
Dot-dashed lines represent an exponential star formation with e- 
folding time r = 12 Gyr, truncated at z =3, with solar metallicity 



values for the local spheroidal luminosity function and by 
Qm and r^A, our data is consistent with no evolution, or 
with modest evolution (at the level of ^ 30%) in terms of a 
decUne in the space density of spheroidals hy z ^ 1. 

• Although we detect little evidence for strong evolu- 
tion in the space density of luminous spheroidals, we find a 
marked deficit in the number of red spheroidals compared to 
predictions where the bulk of star formation was completed 
prior to z :^3. 

• Where redshift data is available it suggests that the 
duty cycle for star- format ion in high-redshift spheroidals is 
indicative of a low level of extended star-formation under- 
lying a passively evolving population, rather than of relic 
star-formation following from a massive burst episode. 

• The apparently mild density evolution and blue 
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APPENDIX A: CALCULATION OF THE 
LUMINOSITY- WEIGHTED MEAN 
METALLICITY OF A SPHEROIDAL GALAXY 
INTERIOR TO AN ISOPHOTAL LIMIT 

Elliptical galaxies are known to contain highly enriched cores 
with strong metallicity gradients. Because of cosmological 
dimming, in the present paper we sample ellipticals at a 
range of rest- frame limiting isophotes, and it is important 
to have at least a qualitative understanding of the effects 
of metallicity gradients on the sampled starlight probed by 
our data. It is shown in Arimoto et al (1997) that the mean 
luminosity-weighted iron abundance of spheroidal systems 
with R^^^ profiles integrated to infinity is similar to the 
abundance measured at the effective radius. It is straightfor- 
ward to generalize the integral in equation (6) of Arimoto 
et al (1997) to calculate {Z{R)), the mean metallicity of 
an R^^^ law spheroid within a given isophotal radius R, as 
a fraction of the metallicity measured within the effective 
radius, Re. For a circularly symmetric galaxy the result- 
ing expression is analytically tractable. The R^^^ law giving 
surface brightness I{R) as a function of radius R is param- 
eterized as follows: 



I{R) = exp < 



{-'[(I) 



1/4 



(Al) 



where Re is the half-light radius, and b = 3.33 In(lO). Metal- 
licity gradients can be parameterized by a power-law index 
c as a function of radius: 



Z{R) = Z(0) 



R_ 

Re 



(A2) 



Essentially all of the ellipticals studied by Arimoto et al. 
(1997) have gradients parameterized by < c < 1.3, with 
metallicities of up to several hundred percent solar in their 
cores decreasing to roughly solar abundance at the effective 
radius. 

Integrating equation (6) of Arimoto et al (1997) to a 
radius R instead of to infinity yields the mean metallicity 
interior to a radius R, expressed in units of the mean metal- 
licity interior to radius Re'. 



(A3) 



{Z{R)) = -e' 
where 

A = [-50406 - 25206^ - 8406^ - 2106^ - 42b^ (A4) 
-76^-6^ + 5040(-l+^^)]. 

[r(8 - 4c) - r(8 - 4c, bR^^^)] (A5) 
B = [-5040(-l + e^^'^^) + 50406i?^/^ (A6) 
+25206^ + 8406^ i?^/^ + 210b^R + 
426' i?'/^ + 76^ i?^/' + b^R^^^] • (A7) 
[r(8-4c) -r(8-4c,6)]. (A8) 

(Note that in the preceeding expression, r(x) is the Euler 
gamma function, and r{x,y) is the incomplete gamma func- 
tion). 

Mean metallicities as a function of radius are shown in 
Figure Al for values of c = 0.4 and c = 1.3, spanning the 
range of metallicity gradients typically seen in luminous el- 
lipticals. As expected, the luminosity- weighted metallicity 
as a function of isophotal radius rises sharply at R <^ Re, 
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Table 1. The HST Sample 



Field 


^814 


Veoe 


H+K' 


FWHM 


a(J2000) 


(5(J2000) 






0029+13 


6300 


3300 


3900 


0.6" 


00:29:06.2 


13:08:12.9 






0144+2 


4200 




2340 


0.8" 


01:44:10.6 


02:17:51.2 






0939+41 


4600 


5400 


3120 


0.9" 


09:39:33.3 


41:32:47.9 






1210+39 


4899 


3999 


3120 


0.6" 


12:10:33.3 


39:29:01.6 






1404+43 


8700 


5800 


3120 


0.7" 


14:04:29.3 


43:19:15.2 






East 


5300 




10985 


0.8" 


12:37:02.0 


62:12:23.4 






HDF 


123600 


109050 


10530 


0.8" 


12:36:47.5 


62:13:04.4 


30 


18 


NEast 


2500 




10595 


0.8" 


12:37:03.7 


62:15:13.0 






NWeast 


2500 




10985 


0.8" 


12:36:49.4 


62:15:53.8 






OutEast 


3000 




10985 


0.8" 


12:37:16.2 


62:11:42.5 






OutWest 


2500 




10270 


0.8" 


12:36:19.3 


62:14:25.5 






SEast 


2500 




10465 


0.8" 


12:36:46.2 


62:10:14.6 






SWest 


2500 




10395 


0.8" 


12:36:32.0 


62:10:55.3 






West 


5300 




11375 


0.8" 


12:36:33.6 


62:13:44.9 






cfrs031 


6700 




3120 


0.8" 


03:02:57.4 


00:06:04.9 


8 


4 


cfrs033 


6700 




3120 


0.7" 


03:02:33.0 


00:05:55.2 


5 


1 


cfrs034 


6400 




3120 


0.8" 


03:02:49.8 


00:13:09.2 


5 


1 


cfrs035 


6400 




3120 


1.0" 


03:02:40.8 


00:12:34.3 


2 


1 


cfrslOl 


6700 




3120 


0.8" 


10:00:23.2 


25:12:45.1 


7 


2 


cfrsl02 


6700 




3120 


0.7" 


10:00:36.5 


25:12:40.1 


9 


6 


cfrsl03 


5302 




3120 


0.6" 


10:00:46.7 


25:12:26.3 


12 


7 


u26xl 


4400 


2800 


4420 


0.8" 


14:15:20.1 


52:02:49.9 






u26x2 


4400 


2800 


3120 


0.7" 


14:15:13.7 


52:01:39.6 






u26x4 


4400 




3120 


0.7" 


14:18:03.2 


52:32:10.4 






u26x5 


4400 




3315 


0.7" 


14:17:56.7 


52:31:00.6 






u26x6 


4400 




3120 


0.5" 


14:17:50.3 


52:29:50.8 






u26x7 


4400 




3120 


0.7" 


14:17:36.9 


52:27:31.2 






u2ay0 


25200 


24399 


4160 


1.4" 


14:17:42.7 


52:28:31.3 






u2iy 


7400 




3120 


0.6" 


14:17:43.1 


52:30:23.3 






u3d3 


4200 


3300 


3120 


1.3" 


20:29:39.6 


52:39:22.3 






ubil 


6300 


3300 


2600 


0.6" 


01:10:00.7 


-02:27:22.2 


5 


1 


ucil 


10800 


4800 


3900 


0.4" 


01:24:41.6 


03:51:24.2 






uehO 


12600 


8700 


3900 


0.5" 


00:53:24.1 


12:34:01.9 


4 


2 


uimO 


11800 




4090 


0.6" 


03:55:31.4 


09:43:31.9 


5 


1 


umd4 


9600 


2400 


3900 


0.7" 


21:51:06.9 


28:59:57.0 






uopO 


4200 


7200 


3120 


0.9" 


07:50:47.9 


14:40:39.1 






uqcOl 


7200 


7800 


3120 


1.0" 


18:07:07.3 


45:44:33.5 






usaO 


6300 


5400 


3120 


0.7" 


17:12:24.0 


33:35:51.2 






usal 


6300 


5400 


3900 


0.5" 


17:12:24.7 


33:36:03.0 








4200 


3300 


3120 


0.9" 


08- 54- 16 9 








ustO 


23100 


16500 


3120 


0.9" 


10:05:45.8 


-07:41:21.0 






ut21 


11999 


3600 


3900 


0.8" 


16:01:13.2 


05:36:03.2 






ux40 


7500 


3300 


3900 


0.6" 


15:19:40.3 


23:52:09.9 






ux41 


6000 


3300 


3120 


0.8" 


15:19:54.6 


23:44:57.8 






uy402 


5400 


1400 


3120 


0.8" 


14:35:16.0 


24:58:59.9 






uzkO 


8700 


8400 


3120 


0.9" 


12:11:12.3 


39:27:04.8 






uzpO 


6300 


3300 


3120 


0.8" 


11:50:28.9 


28:48:35.2 






uzx05 


4700 


2600 


3120 


0.7" 


12:30:52.7 


12:18:52.5 







Table 2. Log of Observations 



Program Date 


Fields Observed 


(Seeing) 


Telescope 


Feb 6-8 1996 


HDF 


1.0" 


UH-2.2m 


Apr 5-8 1996 


HDF 


1.0" 


CFHT 


Apr 17-22 1996 


8 HDF-FF 


0.8" 


UH-2.2m 


Aug 20-24 1997 


21 MDS fields 


0.8" 


UH-2.2m 


Jan 18-20 1998 


23 MDS fields 


0.8" 


UH-2.2m 


Feb 12-13 1998 


11 MDS fileds 


0.7" 


UH-2.2m 



© 0000 RAS, MNRAS 000, 



The Colour Distribution of Faint Field Spheroidals 15 



Table 3. Absolute Numbers 



Data Set 



early 



early 



Total Excess Red Fraction Excess^ 



Observational data (with bootstrapped errors) 



Visually Classified 
Visually Classified 
A/C classified 



Compacts 



266 ± 13 
316 ± 14 
323 ± 13 



26 ±5 
40 ±5 
40 ±6 



POZZETTI ET AL LF^ (a = —0.48) 



40% Solar Qm 
100% Solar Qm 
250% Solar Qm 

40% Solar Qm 
100% Solar Qm 
250% Solar Qm 



40% Solar Qm 
100% Solar Qm 
250% Solar Qm 

40% Solar Qm 
100% Solar Qm 
250% Solar Qm 



= 0.0 Qa = 0.0 


357.17 


72.63 


1 


r = 0.0 Qa = 0.0 


366.57 


135.26 


1 


r = 0.0 Qa = 0.0 


372.85 


198.09 


1 


= 0.3 Qa = 0.7 


433.94 


90.37 


1 


= 0.3 Qa = 0.7 


449.58 


168.27 


1 


= 0.3 Qa = 0.7 


458.74 


247.12 


1 


LF" {a = -1.00) 








= 0.0 Qa = 0.0 


240.29 


28.56 





= 0.0 Qa = 0.0 


259.81 


65.18 





= 0.0 Qa = 0.0 


286.76 


113.26 


1 


= 0.3 Qa = 0.7 


284.40 


35.43 


1 


= 0.3 Qa = 0.7 


311.23 


80.37 


1 


= 0.3 Qa = 0.7 


344.91 


139.39 


1 



1.45 
1.31 
1.47 
1.34 



1 

-'^•"'^1.5 



72 
_.56 
1.78 
1.61 

'^1.64 



nO.95 
^0.86 
Ql.03 
^0.93 



1.08|; 



1.13 
03 



1.12 
1.02 

1 171-23 
-^•-^'1.12 
1 onl.36 
i.C)U^_24 



9 7q3.46 
^- ' ^2.34 

5.206-|4 
7.629^3 



3.48^:30 

6.47i:2J 

9.50^i,r 



1 10^-36 

2.513-10 
4.36t-39 



-'-•'^"1.14 
3.093;83 

5.366-64 



" Single Burst Model with 1.0 Gyr duration and zp = 5. 

^ ^total,model/ ^ early , observed^ bascd on visual classifications with 1(7 upper and lower limits. 

c jyi8i4-HK'>3.0^jylsi4-HK'>3.0 ^^^^^ ^^^^^^ classifications with la upper and lower hmits. 

early ,moaeL ' early , observed ' 



Luminosity Weiglited Metallicity vs. Isophotal Radius 



4 




0.5 1 1.5 2 2.5 3 3.5 4 
R/Re 

Figure Al. The mean luminosity-weighted metallicity as a function of isophotal radius. Metallicity is given as a fraction of the 
metallicity interior to the half-light radius Re. Radii are expressed in units of the half-light radius. The dotted line corresponds to 
c = 1.3; and the solid line to c = 0.4. 
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but for our present purposes the important feature to note 
in Figure 12 is the shallow decline beyond R > 0.5Re in 
the luminosity-weighted mean metallicity, regardless of the 
strength of the power-law slope c parameterizing the metal- 
Ucity gradients. This behaviour was also noted by Arimoto 
et al. 1997 in their simplified calculation. The overall metal- 
licity of the galaxy can be well- characterized by the mean 
metallicity interior to the effective radius (ie. ^ solar in the 
Arimoto et al. sample) over an extraordinarily large range 
in isophotal radius. We can thus make fair comparisons be- 
tween ellipticals over a broad range of redshifts in spite of 
cosmological dimming of the rest-frame isophote, provided 
only the ellipticals are observed with sufficient signal-to- 
noise that their isophotal radii are at least as large as their 
effective radii. This is the case for the vast majority of ellip- 
ticals in the present sample — only at very high redshifts, 
where compact ellipticals have limiting isophotes compara- 
ble to their effective radii, is the mean luminosity-weighted 
metallicity expected to deviate from that at the effective 
radius. 
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